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ABSTRACT 

We have confirmed the existence of the transiting super Neptune extrasolar 
planet HAT-P-llb. On May 1, 2009 UT the transit of HAT-P-llb was detected 
at the University of Arizona's 1.55m Kuiper Telescope with 1.7 millimag rms 
accuracy. We find a central transit time of = 2454952.92534 ± 0.00060 BJD; 
this transit occurred 80 ± 73 seconds sooner than previous measurements (71 
orbits in the past) would have predicted. Hence, our transit timing rules out 
the presence of any large (> 200s) deviations from the ephemeris of Bakos et al. 
(2009). We obtain a slightly more accurate period of P = 4.8878045 ± 0.0000043 
days. We measure a slightly larger planetary radius of Rp = 0.452 ± 0.020-Rj 
(5.07 ± 0.22i?e) compared to Bakos and co-workers' value of 0.422 ± O.OlARj 
(4.73 ± O.I6P0). Our values confirm that HAT-P-llb is very similar to GJ 
436b (the only other known transiting super Neptune) in radius and other bulk 
properties. 

Subject headings: planetary systems, stars: individual: HAT-P-11 



1. Introduction 



The transit of an extrasolar planet across the face of its host star allows direct mea- 
surement of the bulk properties of the planet. In particular, the transit allows accurate 
determination of the planet's radius (see for example Charbonneau et al. 2006 and refer- 
ences within). When these radii are combined with radial velocity (RV) measurements of 
masses then densities of the transiting extrasolar planets can be calculated. Knowledge of the 
heating from the star allows estimates of the temperatures of the irradiated planets. Models 
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of the bulk properties of these planets can be compared to observation (see for example 
Baraffe et al. 2008; Fortney et al. 2007; Burrows et al. 2007; Seager et al. 2007). 

While some 5S0 transiting extrasolar planets are now known, only 3 have masses less 
than 10% of Jupiter (closer to Neptune in mass). The first transiting super Neptune was 
GJ 436b which was discovered around a M2.5 star at 0.028 AU by an RV survey of Butler 
et al. (2004). They found it had a mass of ~21 M®. Follow-up photometric observations 
of GJ 436b then discovered it to be a ~7-8 millimag (~0.7%) transiting planet (Gillon et 
al. 2007). Further follow-up measurements find a radius of ~ 4.2 — 4.9-R© (Torres, Winn & 
Holman 2008; Bean et al. 2008; respectively) and a density of p = 1.69~^Q\2gcm~^ (Torres 

2008) . Through these transit observations and modeling by Baraffe et al. (2008) it has been 
determined that GJ 436b is mainly composed of metals with only a small H/He envelope. 

The second transiting Neptune discovered was HAT-P-llb, in orbit around HAT-P-11 
(2MASS 19505021+4804508) a K4 (V=9.6 mag) metal rich star. This planet was discovered 
by the HATNet array of small 0.11m telescopes on Mt. Hopkins in Arizona (Bakos et al. 

2009) . At just 4.2 millimag (mmag) HAT-P-llb was the smallest planet discovered by the 
transit method. Moreover, Bakos et al. (2009) argue that compared to recent measurements 
of the radius of GJ 436b, HAT-P-llb was at the time the smallest transiting extrasolar 
planet known. However, the CoRoT team has announced the discovery of C0R0T-7b which 
is smaller still with just 1.7 (Rouan et al. 2009). In any case, transiting objects of 
less than 0.1 Mj number no more than three today and start to probe densities and masses 
closer to terrestrial -in contrast to lower density gas giants composed mainly of a large H/He 
envelope. 

In the case of HAT-P-llb, detailed RV measurements by Bakos et al. (2009) show 
a linear drift (0.0297 ± 0.0050 m/s/day) in the RV residual of HAT-P-11. This drift is 
possibly due to the pull of an additional unseen planet in the system (Bakos et al. 2009). In 
addition, Bakos et al. (2009) determine a non-zero (0.198±0.046) eccentricity which might 
be maintained by interactions with another planet. Both observations hint at the presence of 
another outer planet "HAT-P-llc" in the system. Indeed, most systems with super Neptunes 
are multiple planet systems (Bakos et al. 2009 and references within). However, it is worth 
noting that to date no transiting planet is known to be a member of a multiple planet system. 
Hence, detection of multiple transiting planet system would be very interesting. Continued 
RV monitoring of this system may directly detect curvature in the RV residuals due to this 
outer planet. Another way to directly detect the presence of a possible "HAT-P-llc" would 
be a sensitive search for transit timing variations over a series of HAT-P-llb transits. A 



^http:/ /exoplanet.eu 



- 3 - 



search to bound the magnitude of such timing variations motivated this paper. 

2. Observations & Reductions 

Data was taken at the University of Arizona's 61-inch (1.55m) Kuiper telescope on 
Mt. Bigelow, Arizona on 1 May, 2009 UT with the Mont4k CCD, binned 3x3 to 0.43'7pix. 
Observing a single transit of HAT-P-llb at high S/N is fairly challenging, since the transit 
depth is just ~ 4.3 mmag, HAT-P-11 itself is a very bright V = 9.6 magnitude K4 star 
with a nearby faint (likely background) star, the surrounding "Kepler field" is somewhat 
crowded, and all of the potential photometric reference stars are several magnitudes fainter 
than HAT-P-11. Defocusing was not possible, due to the crowded field and the difficulty of 
maintaining a consistent focus offset with this telescope/instrument combination. In order 
to take long enough exposures to sufficiently average out atmospheric scintillation noise, 
while still avoiding saturation of the ccd, we used a medium bandwidth Stromgren b filter 
(AA =18.0 nm). The Mont4k filter holder and filter sensor position were specifically modified 
to accommodate the thicker Stromgren filter for our observations. 

The conditions were photometric with light wind and no moon throughout the observa- 
tional period. In total 448 images were obtained with < 2 pixels of wander, due to excellent 
autoguiding. Integration times of 20 s were used at the start when the target was still at 
relatively high airmass. Just prior to the start of the transit (after the 109*^^ image), the 
exposure times were reset to 17 s, giving a sampling time of 28.3 s for the remainder of 
the observations. The relatively fast overhead time is achieved mainly by binning 3x3 and 
skipping the fiushing of the ccd after each readout and before the subsequent exposure in 
the sequence, but is also a product of the Mont4k's design, which includes two amplifiers 
and pre-amplifiers. 

The images were bias-subtracted, fiat-fielded, and bad pixel-cleaned in the usual man- 
ner. Aperture photometry and sky subtraction was performed for the target star and three 
reference stars using the aperture photometry task PHOT in the IRAF DAOPHOT pack- 
ag^. An aperture radius of 4.3" (10.0 pixels) was adopted, as it produced the smallest scatter 
in the light curve and also eliminated contamination from HAT-P-ll's nearby companion, 
which is 5.4 magnitudes fainter in b and 8.9" to the NNE (at a position angle of 9°). The 
reference stars are 2.3 to 2.65 magnitudes fainter than HAT-P-11, and were chosen to be 



^IRAF is distributed by the National Optical Astronomy Observatories, which are operated by the As- 
sociation of Universities for Research in Astronomy, Inc., under cooperative agreement with the National 
Science Foundation. 
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distributed as uniformly as possible about the target star on the sky (they formed a triangle 
around HAT-P-11 at distances of 89, 208, 345" from HAT-P-11). The reference stars were 
normalized to unity and then weighted according to their average fluxes. 

We applied no sigma chpping rejection to the reference stars or HAT-P-llb -all dat- 
apoints were used in the analysis. The final hght curve for HAT-P-11 was normahzed by 
division of the weighted average of the three reference stars. The residual light curve in Fig. 
1 (bottom left) has a photometric RMS range of 1.7 mmag rms and a time sampling of 29 
seconds on average. This is very typical of the relative photometric precision achieved with 
the Mont4k on the 61-inch (Kuiper) telescope for high S/N images (Randall et al. (2007); 
Dittmann et al. 2009). 

3. Analysis 

3.1. A Search for Transit Timing Vciriations 

The planetary transit hght curves were fit using the method prescribed by Mandel 
and Agol (2002). The transit HAT-P-llb parameters used in the fit were those in Table 1 
measured by Bakos et al. (2009). The correct hnear and quadratic limb darkening parameters 
for our b filter were taken from Claret (2000). In order to detect any transit timing variations 
the only parameter that was allowed to vary in the fit was the central time of of the transit, 
Tc . The time of the center of this transit, (T^), is shown near the bottom of Table 1. We 
note that the purpose of this section of the paper was not to rc-dcrivc all the parameters 
of the transit but to understand if our transit is consistent with the period of Bakos et al. 
(2009). 

To measure we minimized to find a = 2454952.92534 ± 0.00060 BJD value (the 
fiux uncertainty for each datapoint in the fit was calculated from the propagation of the 
photometric errors determined with the PHOT task). The ±0.00060 day la uncertainty was 
estimated by Monte-Carlo simulations of 1000 simulated datasets with the same 1.7 mmag 
rms scatter of the original data (see Fig. 2 left). 

3.2. Determination of a new Planetary Radius Value 

To try and understand if our transit data suggest a different planetary radius for HAT- 
P-11 we repeated the fit in the last section using the method prescribed by Mandel and 
Agol (2002) but this time allowed the planetary radius Rp to vary along with Tc. In Fig. 1 
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(right) we see the result of our fit (sohd red hne) and the residuals of this fit below. 

To measure Rp/ we further minimized the reduced xt to 1.06 with simultaneous fits 
of Rp/R^ = 0.0621 ± 0.0011 and = 2454952.92534. The ±0.0011 la uncertainty in Rp/R^ 
was estimated by Monte-Carlo simulations of 1000 fake datasets with the same 1.7 mmag 
rms scatter as the original data (see Fig. 2 right). 

4. Discussion 

4.1. Is the Timing of the Transits Changing? 

A key goal of this paper is to compare our measured May 1, 2009 UT to that predicted 
from the previously measured values. Projecting the P=4. 8878162 ±0.0000071 day period of 
Bakos et al. (2009) forward from their most accurate = 2454605.89132 ± 0.00032 transit 
suggests that our May 1 transit (n=71 periods later) occurred AT = 80 ± 73 seconds sooner 
(where ctat was calculated by \/(c'"Tci+^'^p+'^rc2)) than our observed Tc2 = 2454952. 92534± 
0.00060 BJD was predicted to be. However, the significance of this disagreement is small. 
Indeed there is ~60% probability that our observations are fully consistent with the timing 
measurements (and uncertainties) of Bakos et al. (2009). Certainly, we can rule out large 
> 200 second timing errors at the ~ 3a level. 

With the addition of our new Tc values to the two previous values we derive a new 
P=4.8878045± 0.0000043 day value for the period of HAT-P-llb (based on a sigma weighted 
average; see Table 1). This new value is slightly shorter than the P=4.8878162 ± 0.0000071 
day period of Bakos et al. (2009). However, it will require future observations to determine 
if this new period will better predict future transit times. It is entirely possible that all of 
our transit timing values are consistent with predictions from Bakos et al.'s ephemeris within 
measurement errors. Moreover, one can calculate that the period estimate between the 
measurements over the first 284 orbits of Bakos et al. (2009) and over the last 71 orbits 
has only changed by 0.177 ± 1.009 seconds. Hence, there is no significant evidence, with the 
data in hand, that Hat-P-llb's period has changed over the last 0.95 year compared to the 
previous 3.79 years. 

In general, it is difficult with a single additional transit to confidently determine if the 
discovery period of Bakos et al. is changing with time. However, our Tc2 datapoint 71 orbits 
later yields a more accurate period of P=4. 8878045 ± 0.0000043 day. This new period with 
its lower error can be used to construct a standard 0-C (observed minus calculated time-of- 
transit) diagram which will start showing "non-linear" trends in the future 0-C residuals if 
the period is truly changing. 
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4.2. What is the Radius of HAT-P-llb? 

Our observations of the depth of the transit finds a deeper transit and a Rp/R^, = 
0.0621 ± 0.0011 compared to 0.0576 ± 0.0009 of Bakos et al. (2009). We find that Rp/R^ 
is 0.0045 ± 0.0014 larger than that of Bakos et al. (2009). Hence, there is a significant 
probability that our transit was deeper than that of Bakos et al. (2009). We derive a 
slightly larger planetary radius of Rp — 0.452 ±0.020 Rj (5.O7±O.22i?0) compared to Bakos 
et al.'s values of 0.422 ± 0.014i?j (4.73 ± O.lGi?®)- Our values suggest that HAT-P-llb is 
very similar to GJ 436b in radius. 

5. Conclusions 

We confirm the existence of the transiting planet HAT-P-llb. Our main conclusions 
from our 1.7 mmag rms (unbinned) transit observations (with the University of Arizona's 
1.55 Kuiper telescope) of the May 1, 2009 UT transit are: 

1. We find a central transit time of = 2454952.92534 ± 0.00060 BJD from a best fit 
to our data. We estimate that the transit occurred 80 ± 73 seconds sooner than previous (71 
orbits in the past) measurements would have predicted (Bakos et al. 2009). Our finding is 
consistent with the ephemeris of Bakos et al. and rules out the presence of any large timing 
variation. 

2. We derive a slightly larger planetary radius of Rp = 0.452 ±0.020 Rj (5.07±0.22i?e) 
compared to Bakos et al.'s values of 0.422 ± 0.014it!j (4.73 ± 0.16i?©). Our values suggest 
that HAT-P-llb is very close to GJ 436b in radius. 

We would hke to thank the anonymous referee for helpful comments leading to a better 
final paper. We would hke to thank Greg Stafford, Gary Rosenbaum, and the Catalina 
Mountain staff for modifying the Mont4k filter box and holder so that we coTild use the 
Stromgren b filter. We would also like to thank the Arizona NASA Space Grant program 
for funding this work. LMC is supported by a NSF Career award and the NASA Origins 
program. 
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Table 1: Parameters of the HAT- P- 11 system 



Parameter Value Reference 

P (days) 4.8878162 ± 0.0000071 Bakos et al. (2009) 

Tc (BJD) 2453217.75466 ± 0.00187 Bakos ct al. (2009) 

Te (BJD) 2454605.89132 ± 0.00032 Bakos ct al. (2009) 

b 0.347i[J:}^^ Bakos et al. (2009) 

i (deg) 88.5 ± 0.6 Bakos et al. (2009) 

(Rq) 0.75 ± 0.02 Bakos et al. (2009) 

Rp/R^ 0.0576 ± 0.0009 Bakos et al. (2009) 

Mp (Mj) 0.081 ± 0.009 Bakos et al. (2009) 

Rp (Rj) 0.422 ± 0.014 Bakos et al. (2009) 

(BJD) 2454952.92534 ± 0.00060 This work 

P (days) 4.8878045^ ± 0.0000043 This work 

Rp/R^ 0.0621 ± 0.0011 This work 

Rp {Rj) 0.452 ± 0.020 This work 



"the new period value was calculated by a sigma weighted least square of all three Tc values 
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Fig. 1.— Left: The timeseries of HAT-P-11 during the Transit of May 1, 2009 UT. We 
show our best fit (reduced xl =1-06) with simultaneous fits of Rp/ = 0.0621 ± 0.0011 and 
Tc = 2454952.92534 ( sohd red curve). The 1.7 mmag rms residuals of the fit are shown 
below. Right: The timeseries of our three calibrator stars (each normalized by the sum of 
the remaining two calibrator stars). The excellent conditions of the night allowed for mmag 
photometry in individual 17 or 20 second exposures even on these fainter reference stars. 
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Measured T_c - Monte Carlo T_c (days) Measured Radii - Monte Carlo Radii (R_p / R*) 



Fig. 2. — Left: One thousand Monte-Carlo (MC) simulations of independent simulated 
datasets following our best-fit transit model each drawn from a population of data with the 
same 1.7 mmag rms as our data in Fig. 1 (left). Based these MC simulations we find with 
our 1.7 mmag rms uncertainty we can constrain to an accuracy of ±0.00060 days (or 
51.84 seconds) at the la level. Right: Here another one thousand MC realizations imply 
the uncertainty in the Rp/R^ ratio to be ±0.0011 at the la level. 



